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Abstract 



We investigate the feasibility of reconstructing the radial intensity profile of ex- 
tended stellar sources by inverting their microlensed light curves. Using a simple, 
linear, limb darkening law as an illustration, we show that the intensity profile can 
be accurately determined, at least over the outer part of the stellar disc, with real- 
istic light curve sampling and photometric errors. The principal requirement is that 
the impact parameter of the lens be less than or equal to the stellar radius. Thus, 
the analysis of microlensing events provides a powerful method for testing stellar 
atmosphere models. 
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1 Introduction 



In recent years a number of authors have investigated the microlensing of ex- 
tended stellar sources. Simmons, Newsam & Willis (1995) (hereafter SNW95) 
have shown that the light curves of extended sources can exhibit a significant 
chromatic dependence, essentially because limb darkening renders the effective 
radius of the star a function of wavelength. Thus, in addition to improving 
constraints on the lens parameters, modelling the microlensing of extended 
sources provides a powerful tool for gravitational imaging stellar surfaces. 

In this contribution we describe how microlensing light curves of extended 
sources may be used as a test of stellar atmosphere models. We generate 
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artificial light curves, assuming a particular model atmosphere, and use the 
Backus- Gilbert numerical inversion method to estimate the radial stellar in- 
tensity profile from the observed light curves. 



2 Inverting Microlensing Light Curves 

The (time dependent) integrated flux, F(t), from an extended stellar source 
of radius, R, lensed by a point lens is given by (c.f. SNW95) 

R 2tt 

F{t)= J J I(s,0)A{d)sdsdO (1) 

s=0 6»=0 

where d is the projected distance from the lens to the element of the stellar 
surface and the amplification function, A(d), takes its usual analytic form for 
a point source. (Note that d is a function of s, 9 and t). 

If we assume that the projected stellar surface displays circular symmetry - 
i.e. I(s, 9) = I(s) - we may write the above equation in the form 

i 

F(t) = J I(r)K(r,t)ds (2) 

r=0 

where r = s/R and the kernel function, K(r, t), is obtained by integrating over 
9. A solution to this integral equation for I(r) can be obtained by applying 
the Backus-Gilbert inversion procedure. This method takes account of the 
smoothing effect of the kernel function and reconstructs a regularised estimator 
of I{r) which optimises the trade-off between the bias and variance of the 
estimator. For details of the method see Coleman, Gray & Simmons (1998) 
and Coleman (1998). 



3 Testing Stellar Atmosphere Models 

In order to test the feasibility of reconstructing stellar surface profiles, we 
assumed a simple, linear, limb darkening law with coefficient, u. We generated 
microlensed light curves of typically 100 data points, with an even sampling 
rate, and Gaussian noise added to the photometry at a level of typically 2% of 
the baseline flux. In most cases we took the impact parameter and the einstein 
radius equal to the stellar radius. We considered two model atmospheres: 
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(1) T cS = 4000/T; logg = 2.0; u B = 0.92, u v = 0.82 



(2) T cS = 6730fT; logg = 4.5; u B = 0.76, u v = 0.63 

where the Johnson B and V band linear limb darkening coefficients are from 
Claret & Gimenez (1990). We carried out inversions of the B and V band 
intensity profiles for these models, and investigated the effect of changing the 
impact parameter, stellar radius, light curve sampling rate and photometric 
accuracy. 



4 Results 



Figure 1 illustrates the reconstructed V band profile for model (2), for the 
case of the impact parameter and stellar radius equal to the einstein radius 
of the lens. The errors on the recovered solution are determined from the 
covariance matrix of the Backus-Gilbert estimator. We can see that the recon- 
structed profile is significantly biased for r<, 0.6, due to the smoothing effect 
of the kernel function, but the true profile is well recovered over the interval 
0.6 <r <, 0.9. We now briefly summarise the results of varying the stellar and 
lens parameters. 

• The inversions generally recover the true profiles well over the interval 
0.6 < 0.9, for a wide range of T e $ , log g and u. 

• Increasing the einstein radius improves significantly the inversions for r ■> 0.6, 
but no improvement is seen for r £ 0.6. 

• Reducing the impact parameter (i.e. a transit event) significantly improves 
the accuracy of the reconstruction for r > 0.6 and reduces the bias for r <, 0.6. 
For impact parameters greater than the stellar radius, however, the recon- 
struction deteriorates rapidly for all r. 

• Even with photometric errors of 10% a reasonable recovery of J(r) is still 
obtained over the interval 0.6 <r £0.9; on the other hand, reducing the 
errors to only 0.2% does not improve the recovery for r £ 0.6, however. This 
is because the bias is primarily due to the ill-posedness of the kernel over 
this range, and not due to the photometric errors. 

• The reconstructions become unacceptably noisy when the number of light 
curve data points is reduced to n — 10, but there is little further improve- 
ment in accuracy above n = 50. 
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Fig. 1. Recovered radial V band intensity profile from a light curve of 100 points. 
The true profile, for uy = 0.63, is shown by the dotted line 

5 Conclusions 



Our results indicate that - with realistic light curve sampling and photometric 
errors - one can accurately reconstruct, at least in part, the multicolour ra- 
dial intensity profiles of extended stellar sources from their microlensed light 
curves, provided that the impact parameter of the lens is comparable to the 
stellar radius. The smoothing properties of the kernel function result in a bi- 
ased solution for r < 0.6, unless the lensing event is a transit with small impact 
parameter. Nevertheless, the accurate recovery over the interval 0.6 <,r £0.9 
is a robust result over a wide range of stellar temperatures and limb darken- 
ing coefficients. Despite the narrow width of this 'good fit' annulus, it is still 
adequate to usefully discriminate between different model atmospheres - e.g. 
two models with the same temperature but with log g = 2.0 and log g = 4.0. 
Thus, we conclude that broad band microlensed photometric light curves are 
a powerful tool for investigating extended stellar sources and testing model 
stellar atmospheres, and form a useful adjoint to spectroscopic and polarimet- 
ric microlensing signatures. We are currently investigating the application of 
inversion techniques to more realistic model atmospheres and stellar intensity 
profiles. 
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